UNCLASSIFIED 


_ AD  NUMBER _ 

AD843091 

LIMITATION  CHANGES 
TO: 

Approved  for  public  release;  distribution  is 
unlimited. 


FROM: 

Distribution  authorized  to  U.S.  Gov't,  agencies 
and  their  contractors ; 

Administrative/Operational  Use;  JUN  1968.  Other 
requests  shall  be  referred  to  Naval  Ship 
Research  and  Development  Center,  Washington, 

DC. 


_ AUTHORITY 

NSRDC  ltr  12  Oct  1972 


THIS  PAGE  IS  UNCLASSIFIED 


AD843091 


FURTHER  CALCULATIONS 
OF  THE  FLUTTER  SPEED  OF  A  FULLY 
SUBMERGED  SUBCAVITATING 

HYDROFOIL 

by 

Wan-Hwo  Chu  and  H.  Normoi.  Abramson 


Technical  Report 
Contract  No.  N00014-68-C-0259 
SwRI  Project  02-2311 


Prepared  for 

The  U.S.  Navy 

Naval  Ship  Research  and  Development  Center 
Washington,  D.C.  20007 


June  1968 


This  research  mm  earned  out  under  the  Naval  Ship  System,  Command 
Hydrofoil  Exploratory  Development  Program  SF  013  03  01,  administered 
by  the  Naval  Ship  Research  and  Development  Center.  Prepared  under  the 
Office  of  Naval  Research  Contract  N00014-WC-03S9 


This  document  is  subject  to  special  export  controls  and  each  transmittal  to 
foreign  governments  or  foreign  nationals  may  br  made  only  with  prior 
Mprovel  of  the  Head,  Hydromechanics  Laboratory.  Naval  Ship  Research 
and  Development  Center.  Washington.  D.C.  30007. 


SOUTHWEST 

SAN  ANTONIO 


RESEARCH  INSTITUTE 

HOUSTON 


aSttHHMBl 


SOUTHWEST  RESEARCH  INSTITUTE 
8500  Culobro  Road,  San  Antonio,  Texas  78228 


FURTHER  CALCULATIONS 
OF  THE  FLUTTER  SPEED  OF  A  FULLY 
SUBMERGED  SUBCAVITATING 

HYDROFOIL 

by 

Wen-Hwo  Chu  and  H.  Norman  Abramson 


Technical  Report 
Contract  No.  N00014-68-C-0259 
SwRI  Project  02*2311 


Prepared  for 


The  U.S.  Navy 

Naval  Ship  Research  and  Development  Center 
Washington,  D.C.  20007 


D  D  C 

tm  npfj? 

.MOV  18 as  |i 

IJEtSKu  U  l±® 

D 


June  1968 


This  research  mm  carried  out  under  the  Neva!  Ship  Syttems  Command 
Hydrofoil  Exploratory  Development  Program  SF  013  0301,  administered 
by  the  Naval  Ship  Reeearch  end  Development  Center.  Prepared  under  the 
Office  of  Naval  Research  Contract  N0001 4-68 -C -0259. 


This  document  is  subject  to  special  export  controls  end  each  transmittal  to 
foreign  governments  or  foreign  nationals  may  be  made  only  with  prior 
approval  of  the  Head,  Hydromechanics  Laboratory,  Naval  Ship  Research 
and  Development  Canter,  Wellington,  D.C.  20007. 


Approved:  y 

— . 

H.  Norman  Abramson,  Director 
Department  of  Mechanical  Sciences 


ABSTRACT 


Some  new  calculation*  of  the  flutter  speed  of  the 
SwRI  fully  submerged  subcavitating  hydrofoil 
model  are  presented.  Variations  in  lift  curve 
slope  and  in  center  of  pressure  location  are  found 
to  have  a  most  profound  influence  on  both  flutter 
speed  and  frequency.  When  variations  in  these 
parameters  are  combined  with  a  relaxation  of  the 
Kutta  condition  (proposed  previously),  excellent 
agreement  with  the  measured  flutter  speed  is 
obtained. 
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INTRODUCTION 


The  catastrophic  structural  failure  of  the  SwRI  flutter  model  as 
reported  by  Abramson  and  Ransleben  (Ref.  1)  has  been  the  focal  point  of  much 
discussion  and  the  subject  of  a  variety  of  calculations  and  analyses.  The  con* 
ditions  at  which  catastrophic  failure  of  this  model  occurred  were  a  velocity 
of  48. 1  knots  and  a  frequency  of  17.  5  cps  (corresponding  to  a  reduced  veloc¬ 
ity  of  1.  48).  However,  in  making  any  comparison  or  correlations  with  these 
data,  it  should  be  quite  clearly  noted  that  uncertainties  exist  even  here 
(Refs.  2,  3),  so  that  exact  agreement  should  be  viewed  with  as  much  skepticism 
as  would  wide  disagreement. 

Initial  attempts  to  calculate  the  flutter  velocity  of  this  model  by  two* 
dimensional  and  Reissner-Stevens  theories  failed  to  yield  reasonable  results 
(Refs.  1,  4),  although  the  model  had  an  aspect  ratio  of  five.  This  difficulty  was 
not  unanticipated,  however,  as  large  differences  between  experimental  data  and 
computed  results  for  various  flutter  studies  at  low  mass  density  ratios  had 
previously  been  encountered  and  had  been  the  subject  of  much  controversy 
(Refs.  2,  5).  It  therefore  has  now  become  rather  widely  accepted  that,  for 
the  range  of  parameters  pertinent  to  hydrofoil  applications  the  more*or*less 
classical  formulations  of  the  unsteady  hydrodynamic  loads  are  not  adequate. 
However,  even  the  evidence  available  here  is  somewhat  clouded:  for  example, 
measured  data  on  oscillatory  lift  and  moment  distributions  (Refs.  6,  7)  employed 
directly  in  a  flutter  analysis  of  the  SwRI  flutter  model  yielded  a  substantially 
overconservative  prediction  (Refs.  1,  4);  unfortunately,  the  accuracy  of  this 
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measured  data  is  also  open  to  some  question  (Ref.  3),  and  the  analytical 
techniques  employed  in  the  analysis  may  not  have  been  adequate. 

A  somewhat  different  approach  was  also  taken  quite  early  in  all  of 
these  studies  by  considering  some  relaxation  of  the  Kutta  condition  (Refs.  8,9); 
but*  again,  reasonable  flutter  predictions  were  not  readily  forthcoming  (Ref. 4) 
and  could  only  be  obtained  from  this  semiempirical  procedure  by  using  rather 
large  arbitrary  shifts  in  the  phase  of  the  circulation  function. 

Several  years  ago,  Yates  began  the  development  of  a  modified  strip* 
analysis  method  employing  arbitrary  spanwise  distributions  of  lift  curve  slope 
and  center  of  pressure  location  (Refs.  10-14).  This  method  can  be  based  on 
the  idea  of  employing  measured  data  for  these  parameters  or,  alternatively, 
employing  values  calculated  from  lifting  surface  theory.  Yates  has  applied 
this  basic  method  to  the  SwRI  flutter  model  (Ref.  15),  obtaining  a  flutter 
velocity  of  approximately  38.  5  knots,  with  a  structural  damping  coefficient 
of  g  =  0.  02  cps  and  a  flutter  frequency  of  21, 3  cps.  Some  further  refinements 
and  modifications  appear  to  have  yielded  somewhat  better  agreement,  but 
these  have  not  yet  been  reported  in  detail. 

Yet  another  approach  has  relied  upon  the  rather  full  application  of 
lifting  surface  theories  (Refs.  16-18).  The  method  proposed  by  Rowe  (Ref.  18) 
provides  for  some  variation  in  satisfaction  of  the  Kutta  condition  and  inherently 
generates  data  on  lift  curve  slope  and  center  pressure  location. 


It  it  the  purpose  of  the  present  report  to  present  yet  another  series 
of  results.  These  are  all  based  essentially  upon  strip  theory,  but  Involving 
various  modifications  to  account  for  Kutta  condition,  lift  curve  slope, 
center  of  pressure  location,  etc. 
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ANALYSIS 


Basic  Equations 

For  combined  pitching  and  heaving  oscillations  of  a  hydrofoil, 
represented  as  a  uniform  unswept  beam,  the  equations  of  motion  may  be 
put  into  the  nondimen sional  form  (Ref.  19) 

u2 

-HYflhJ  -  Hyahxaa0  +  -f  0  +  **>  *»o  *  AShh0  +  A£a°0  <U> 

J- 

T‘1fahxah0  •  l‘r!l'iaO  +  I'l  J  ll  +  *r&a  *  AJhhS  +  A!>0  <lb> 

where 

r 1  -1  -1 

7hS  J  fh(y**  dy*  ’  y\  S  J  dy*  *  Yoh  *  /  fa(y*>  fh(y*>  dy* 

(2a,  b,  c) 

f&,  fjj  are  the  pitching  and  heaving  mode  shapes,  respectively,  i.  e. , 


a  *  aofafy*> 

positive  nose-up 

(3a) 

£sh0fh(y*> 

positive  downward 

(3b) 

y*  s  X.  |  = 
7  sb 

aspect  ratio,  b  =  semichord 

(4a) 

„  M 

p  * - r  . 

2  la 

ra  =  2 

(4b,  c) 

irpb2 

Mb2 

pZ.4 

(4d.e) 

5 


The  (nondimensional)  generalized  loadings  are 

,1 

Ahh*/  Lhfh<y*>dy* 

0 

,1 

Aho  *  /  LS£a<y*>  fh<y*>  dy* 
j 

Kh*  I  Mhfh(y*)  fo(y*)  dy* 

0 

Aja  -/  M*fJ(y*)dy* 

0 


(5a) 

(5b) 

(5c) 

(5d) 


The  force  L*  is  positive  downward;  the  moment  M*  is  positive  nose-up. 

The  determinant  of  the  coefficients  must  vanish  at  the  critical  (flutter) 
speed;  in  this  case. 

agz2  +  ajz  +  a0  =  0  (6) 

where 


z  = 


-f  (1  +  iS) 

w4 


and 

2  2  2  2  2 
*2*1*  rurw?aYh 

ai s  +  KJ  -  yl»rl + 


(7) 

(8a) 

(8b) 


a0  s  +  Ahh>  krhl  +  ASa>  -  fcVah  +  Aha>  ^Vah  +  AJh>  <8c> 
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The  standard  procedure  employed  for  determining  the  flutter  speed  from 

Eq.  (6)  is  to  assume  a  value  for  the  reduced  frequency  k  =  ~and  then  to 
c^- 

cal  culate  z  and  ~  .  The  value  of  g  may  then  be  found  and  a  plot  marie  of 
V  vs.  g;  the  value  of  V  at  which  g  *  0  then  gives  the  flutter  velocity. 

Case  1.  Reissner-Stevens  Method  (Ref,  201 

Following  the  convention  of  Reference  19,  the  expressions  for  force 
and  moment  corresponding  to  pitch  and  heave  are,  respectively! 


r  * 
Lu  = 


L*  s 
a 


h*  . 

2 

r 

irpw^b^ 

"k2 

l 

a 

2 

ii 

irpw^b2 

«•  mmmm  « 

k2' 

L2 

£  +  ik  [C(k)  +  «rh]}  *LjeiOLh 


ji(ik+k2a)+  |\+ik(±-a)j  [C(k) 


K 


w*  _  h 
Mh  " 


M' 

Q 

q  , 


' ${*!■• (■  * i) 11 'clkl •  •i'1} 


*  [l +ik  a)J  [C(k)  +  <ro]  J  *  Mje^o 


(9d) 


The  procedure  required  to  calculate  the  finite  span  correction  factors 
and  <ra  are  given  in  Reference  20  (for  two-dimensional  theory,  s  =  0). 

The  lift  and  moment  calculated  for  the  SwRI  flutter  model  are  found 
to  be  in  agreement  with  those  obtained  previously  with  a  different  computer 


program  (Ref.  1).  The  V-g  curve  is  given  in  Figure  1;  the  flutter  speed  lies 
between  two  previously  calculated  values,  all  of  which  are  based  upon  certain 
data  read  from  curves  given  in  Reference  20.  so  that  the  basic  computer 
program  to  be  employed  in  the  calculation  of  the  curves  to  follow  is  confirmed 
to  be  correct. 

The  measured  mode  shapes  and  characteristics  of  the  SwRI  flutter 

model*  as  employed  in  these  calculations,  are  given  in  Table  1. 

Case  2.  Calculation  Using  Experimental  Data  for  Steady-State  Center  of 
Pressure  Location 

The  expressions  for  lift  and  moment  are  now  written  as 


Lh-  --^{-1^  +  Ik  (c<k) +»!,]}  (10.) 

lJ.  -^{i(lktk*.)+  [l  +  ik(j-  »)]  [C(k)+  <r0]  j.  (10b) 

mJ  *  '  ^  {  2  ’  (*  '  xcp> ik  [C<k>  +  'hi}  (10c) 

Mi  =  -  ^{i  *k  (z  '  m)  -  i  <k2  +  -  <»  -  *cp>  [l+‘k(M] 

•[C(k)+«ra]}  (10d) 


where  xcp  is  the  location  of  the  center  of  pressure,  dependent  upon  the  span- 
wise  positions,  or  obtained  from  the  experimental  data  of  Reference  6.  In 
the  moment  expressions,  Eqs.  (10c)  and  (lOd),  the  term  a  -  xcp  has  been 
used  in  place  of  the  term  -j  +  a  as  the  lift  force  is  to  be  taken  as  acting  at  xc^ 
rather  than  at  the  quarter-chord  point  jbV 


The  V-g  curve  is  given  in  Figure  2,  indicating  a  flutter  speed  of 


approximately  67  knots  (compared  with  the  experimental  value  of  48  knots). 

Case  3.  Calculation  Using  Experimental  Data  for  Oscillating  Center 
of  Pressure  Location 

The  forces  will  be  considered  to  be  the  same  as  those  of  the  classical 
theory  of  Case  1.  but  the  moments  will  be  taken  as 

Mh  *  Mf>iaMh  .  M*  co.  (aMh  -  aLh)  «ioU>  +  (1  L)t 

M*  *  .  M*  co.  <aMo  -  «L(1)  e^Lh  +  x*  Lje^Lh  (ub) 

rQ 

fc  * 

The  unsteady  center  of  pressure  locations  are  taken  as  x__  ,  x_ _  in  these 

c"h  c"a 

expressions  and  defined  (Ref.  6)  so  that 

LQxcpQ  *  Ma  «•  <aMh  -  “Lh>  (12») 

^p,,  *  Mh  co>  (“Mo  ‘  “La1  <12b> 

where 


=  2x 


1 

cPh'  2 


d 


(13a) 


1 

2 


d 


(13b) 


and  xc«  ,  xc_  are  the  oscillatory  center  of  pressure  locations  corresponding 
mi 

to  heave  and  pitch  motions,  respectively,  in  fractions  of  chord  measured 
from  the  leading  edge;  d  is  the  distance  between  the  moment  axis  and  the 
quarter-chord  point.  For  the  SwRI  flutter  model 


TMJ.  m;, 

aMh’  aMa'  Lh-  Li-  °Lh»  aLa  are  given  by  Eqs.  (9a -d)  with 


“h  *  'o  *  °- 


X  +  d  =  0.  3  1 
2 
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(14) 


Values  of  x’’'  ,  x*  are  given  in  Table  2. 

'-rh  '-Pa 

The  V-g  curve  is  given  in  Figure  3,  indicating  a  flutter  speed  of 
approximately  70  knots. 

Case  4.  Calculation  Using  Modified  Values  of  Lift  Curve  Slope  and 
Center  of  Pressure 

It  is  well  known,  in  steady  flow,  that  the  actual  lift  curve  slope  of  a 
foil  is  somewhat  less  than  the  theoretical  value  of  2tt,  thus  having  the  effect 
of  reducing  the  value  of  the  circulation.  Denoting  the  ratio  of  actual  to 
theoretical  lift  curve  slopes  by  (3,  the  Theodorsen  function  may  be  multiplied 
by  this  ratio  (no  additional  corrections  are  to  be  applied  to  the  Reissner- 
Stevens  three-dimensional  factor).  Then,  considering  the  lift  force  to  act 
at  the  center  of  pressure,  the  expressions  for  lift  and  moment  become 


w 

1 

II 

* 

^.id+ik  [pC(k)  +  trh]| 

(15a) 

L*=  -  J-i 

Q  k2  < 

±(ik  +  k2a)+  [l  +ik(±-  a)]  [pC(k)  +  cr0] 

[  (15b) 

mr*  -  2 

Mh  =  -  “7  ■ 
k2 

(a  -  xcp)  ik  [pC(k)+(rh]J 

(15c) 

m:  ■  •  # 

[Mi-*)  -ik2(i+l2)-(a-xcP)  [> 

*»(H 

•  [pc<k)  +  <r0]j 


<15d) 


Various  values  of  C/q  and  xCp  are  given  in  Table  3.  The  data  in  Table  3a  are 
derived  from  experiments  (Ref.  6)while  those  in  Table  3b  are  derived  from 
lifting  surface  theory  (Ref.  15). 


V-g  curves  calculated  by  using  these  two  sets  of  Cja  and  xCp  values 
are  given  in  Figures  4a  and  4b.  The  first,  employing  measured  foil  data,  gives 
a  flutter  speed  of  approximately  100  knots  while  the  second,  employing 
calculated  foil  data,  gives  a  flutter  velocity  of  approximately  76  knots. 
Comparing  these  results  directly  with  those  of  Case  2,  it  would  mean  that 
the  effect  of  modifying  the  lift  curve  slope  in  this  manner  is  to  raise  the 
flutter  velocity. 

Case  5.  Calculation  Using  Yates1  Modified  Strip  Theory  (Ref.  10) 

According  to  Yates  (Ref.  10),  the  constant  factor  i-  -  a  in  the  downwash 
terms  of  a  conventional  strip  theory  analysis  could  be  replaced  by 

X  =  P  +  xcp  -  a  (16) 

which  varies  along  the  span  of  the  foil.  Since  this  modification  includes  three- 
dimensional  effects,  <ra  and  0^  are  to  be  taken  as  zero.  The  expressions  for 
lift  and  moment  then  become 


T  *  -  2 

{-l£  +  ikC(k,} 

(17a) 

L*=  --SJ 
°  k2  1 

'  i(ik+k2a)+  (1  +  ik\)  C(k)j 

(17b) 

<=-  \ 
k  4 

(*?-(£♦•)  ikC(k>} 

(17c) 

(17d) 


1 1 


M°=  •^{lik'-7k2(iti‘2)  '  («  +  £)“  ♦,kMc<k>} 


Two  subcases  art*  now  lonsidered.  Case  5a  is  based  on  the  use  of 
the  values  derived  by  Yates  (Ref.  15)  and  given  in  Table  3b.  The  V-g  curve 
corresponding  to  this  calculation  is  shown  in  Figure  5a  and  gives  a  flutter 
velocity  of  approximately  30  knots  with  g  =  0  and  32  knots  with  g  =  0.  02;  this 
compares  with  a  value  of  38.  5  knots  with  g  =  0.  02  as  computed  by  Yates 
himself  (Ref.  15).  Note,  however,  that  the  slope,  of  the  critical  V-g  curve 
for  small  values  of  g  is  so  flat  that  only  very  slight  changes  in  the  computed 
points  would  be  necessary  to  alter  the  flutter  velocity  by  a  very  large  amount. 
Further,  Yates  employed  six  calculated  modes  (two  bending  and  four  torsion) 
in  his  analysis,  while  the  present  calculations  are  based  on  two  measured 
mode  shapes  (one  bending  and  one  torsion)*.  In  any  event,  it  would  appear 
that  the  modified  strip  theory  proposed  by  Yates  (and  as  employed  in  this 
report)  is  overconservative  by  a  significant  degree. 

Case  5b  employs  the  measured  steady  state  values  shown  in  Table  3a. 
The  corresponding  V-g  curve  is  shown  in  Figure  5b,  giving  a  flutter  speed  of 
approximately  26  knots. 

Case  6.  Calculations  Based  on  Appliedatlon  of  the  Generalized 
Kutta  Condition 

As  proposed  a  number  of  years  ago  (Refs.  4,  8,  9,  21),  a  generalization 
of  the  Kutta  condition  can  be  formulated  in  terms  of  a  factor  applied  to  the 

*  Some  rather  simple  considerations  show  that  our  use  of  the  two  measured 
modes  should  give  reasonably  accurate  results.  Any  discrepancies  arising 
between  using  calculated  and  measured  modes  is  the  result  of  differences 
in  the  dominant  terms  rather  than  in  an  insufficient  number  of  modes. 


trailing  edge  tangential  velocity.  Defining  Kj  =  1  -  rje1^,  and  replacing  the 
constant  factor  ^  t  a  by  a  -  xCp  as  was  done  in  Case  2,  the  lift  and  moment 
become: 


*  2 
Lh  =  -  -4  < 

k2 

(C(k)  +  «rh]  J 

(18a) 

p 

i  (ik  +  k2a)  + 

.  ^ 

[,  +  ,k(H] 

nei4>0  (C(k)  +  <rQ]J. 

(18b) 

2 

Mh  =  -  ”7  ' 
k2 

(*2“  "  (a  "  X°P 

)  ikne^O  [c(k)  + 

-1  - 

(18c) 

(*♦■*) 

(a  *  xcp) 


[\+ik(±-a)j  r?ei<l>0[c(k)+  aa]J  (18d) 


We  shall  assume  that 


n  = 


P 

cos  <j>0 


P  * 


(19a) 

(19b) 


and  shall  select  specific  values  of  $q. 

Case  6a  is  based  on  the  use  of  measured  data  for  Cj^  and  xCp  (Table  3a), 
with  V-g  curves  corresponding  to  several  values  of  shown  in  Figure  6a. 

The  value  of  4>q  *  0  gives  a  flutter  speed  of  approximately  42  knots,  compared 
with  the  experimental  value  of  48.  1  knots;  a  calculation  using  the  value  of 
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4>0  =  10*  gives  almost  exactly  48  knots.  It  is  interesting  to  note  that  the 
Kutta  condition  modification  alone  (e.  g. ,  without  any  simultaneous  correction 
on  CiQ  and  xCp)  gives  computed  results  that  agree  with  the  actual  hydrofoil 
flutter  velocity  of  48  knots  only  for  the  somewhat  unrealistic  value  of  $q  =  -30* 
(Ref.  4).  Even  for  4>o  =  0,  however,  these  new  results  are  quite  good,  especially 
if  one  considers  that  we  have  spoken  only  of  the  flutter  speed  corresponding  to 
g  =  0;  in  fact,  based  on  zero  forward  speed  tests,  the  damping  of  the  SwRI 
flutter  model  was  estimated  to  be  g^  =  0.  016  and  gQ  =  0.  042  in  air,  and  is 
roughly  estimated  to  be  g^  =  0.  070  and  gQ  =  0.  047  in  water,  t  Calculated 
data  corresponding  to  the  critical  frequency  are  also  shown  in  Figure  6a,  giving 
a  value  of  —  *  0.  82  and  j-  *  1. 34  (for  g  -  0)  which  are  in  good  agreement  with 

<*>Q  K 

the  measured  values. 

Case  6b  is  based  on  the  use  of  Yates'  calculated  data  for  C|&  and  xCp 
(Table  3b),  with  V>g  curves  corresponding  to  several  values  of  $q  shown  in 
Figure  6b,  along  with  data  for  critical  frequency.  The  value  of  $o  “  0 
(for  g  =  0)  gives  a  flutter  velocity  of  45  knots  and  -  0.  82  (  =  1.43J  . 

These  are  again  in  very  good  agreement  with  the  measured  values. 

Finally,  Case  6c  differs  from  Case  6a  only  in  that  the  center  of 
pressure  is  taken  as  xCp  =  -  j  ,  with  the  results  being  shown  in  Figure  6c. 

Again,  with  =  0  and  g  =  0,  the  flutter  speed  is  45  knots  and  —  =  0.  82 


tData  provided  by  Mr.  Guido  E.  Ransleben,  Jr. 
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DISCUSSION  AND  CONCLUSIONS 

A  summary  of  the  results  obtained  from  the  various  cases  studied  is 
given  in  Table  4.  Of  these,  only  Case  6  gives  results  that  are  closely 
comparable  with  the  measured  values  of  flutter  speed  and  frequency. 

There  can  be  no  question  that  variations  in  lift  curve  slope  and  in 
center  of  pressure  location,  as  suggested  by  Yates  (Ref.  15),  have  a  most 
profound  influence  on  the  flutter  speed  and  frequency  of  a  foil  having 
parameters  similar  to  those  of  the  SwRI  flutter  model  (Ref.  1).  While  the 
calculated  results  presented  in  this  report  are  not  in  exceptionally  close 
agreement  with  those  of  Yates,  it  appears  this  may  be  the  consequence, 
in  part,  of  sensitivity  of  the  computations  and  also  of  the  use  of  measured 
modes  in  the  present  calculations  and  of  calculated  modes  by  Yates.  Both 
sets  of  results,  however,  load  one  to  conservative  predictions,  and  this  in 
itself  is  a  significant  advance  over  practically  all  previous  results.  The 
difference  between  using  the  measured  values  of  lift  curve  slope  and  the 
center  of  pressure  location  and  that  of  using  those  calculated  from  lifting 
surface  theory  is  apparently  not  large.  Yates'  method  certainly  deserves 
further  study  and  application. 

The  semiempirical  method  based  on  generalization  of  the  Kutta 

conditions  (Refs.  8,  9,  21 ),  and  used  here  in  conjunction  with  lift  curve  slope 

and  center  of  pressure  modifications,  appears  to  give  excellent  results  (Case  6). 

As  with  the  Yates'  method,  the  differences  arising  from  various  Ci  and  xr_ 

a  '-r 

values  are  small.  An  essential  point  is  again  that  the  flutter  predictions  are 
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conservative,  as  well  as  being  in  close  agreement  with  the  measured  flutter 
speed  and  frequency.  It  would  mean  that  the  applications  of  the  procedures 
employed  in  Case  6  to  other  hydrofoil  flutter  studies  would  be  quite  useful 
and  revealing,  if  done  prudently. 

The  role  of  structural  damping  could  be  significant;  in  Case  6,  the 
inclusion  of  small  positive  damping  would  serve  to  improve  even  further  the 
agreement  between  calculated  and  measured  flutter  speed. 
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TABLE  la.  SwRI  FLUTTER  MODEL-MEASURED 
MODE  SHAPES  (IN  AIR) 


y/L 

fh 

*a 

Root 

0 

0 

0 

0.1 

0.025 

0.213 

0.2 

0.075 

0.403 

0.3 

0.158 

0.  570 

0.4 

0.260 

0.710 

0.5 

0.370 

0.825 

0.6 

0.488 

0.915 

0.7 

0.610 

0.966 

0.8 

0.738 

0.990 

0.9 

0.870 

0.999 

Tip 

1.0 

1.000 

1.000 

i 

I 


TABLE  lb.  SwRI  FLUTTER  MODEL 
CHARACTERISTICS  (Ref.  1) 


_ _ Model  Parameter 

Aspect  ratio 
Semichord 

Elastic  axis  location,  a 

Center  of  gravity  location,  xa 

Radius  of  gyration,  r£ 

Bending  stiffness*,  El 

Torsional  stiffness*.  CJ 

Frequency  ratio,  u./u 

n  a 

Torsional  frequency,  ua 
Total  weight  (wing  only) 


Measured  Value 
5.00 
0.  50  ft 
-0.50 

(from  a.  e.)  0.524 

0.512 

3.40  X106  lb-in* 
0.973  X106  lb -in* 
0.490 
20. 5  cps 
121.2  lb 


!?/•  ?a1!Uh!fAen  W4anantWh  arC  ab°Ut  18<?  Cps  and  U*2  CP“-  The  calculated 
about  O.W.  However,  the  measured  frequencies  are  possibly  more 

ccurate  than  the  calculated  values  based  on  measured  El  and  GJ. 


.  * 
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TABLE  2.  VALUES  OF  OSCILLATORY  CENTER  OF  PRESSURE 
LOCATIONS  IN  SEMICHORDS  (Ref.  6) 


1/k 

k 

2xcph 

2iCP„ 

xcph 

x* 

XCpQ 

0.4 

2.5 

0.85 

0.92 

0.54 

0.61 

0.8 

1.25 

0.630 

0.76 

0.32 

0.45 

1.0 

1.00 

0.  560 

0.70 

0.25 

0.39 

1.5 

0.667 

0.480 

0.58 

0.  17 

0.27 

1.8 

0.  556 

0.48 

0.53 

0.  17 

0.22 

2.0 

0.5 

0.48 

0.  50 

0.  17 

0.19 

TABLE  3a.  EXPERIMENTAL  VALUES  FOR  STEADY 
STATE  LIFT  CURVE  SLOPE  C*  AND  CENTER  OF 
PRESSURE,  xcp  (ifef.  6) 


y* 

Cja  (rad.) 

xCp  (semichords) 

0 

4.01 

-0. 522 

0.1 

4.01 

-0.522 

0.2 

4.01 

-0.536 

0.3 

4.01 

-0.  540 

0.4 

4.01 

-0.544 

0.5 

4.01 

-0. 552 

0.6 

4.01 

-0.564 

0.7 

3.724 

-0. 580 

0.8 

3.438 

-0.596 

0.9 

2.  58 

-0.632 

TABLE  3b.  THEORETICAL  VALUES  FOR  STEADY 
STATE  LIFT  CURVE  SLOPE  Cfa  AND  CENTER  OF 
PRESSURE  xcp  FOR  THE  SwRI  FLUTTER  MODEL 


y* 

C ia  (rad.) 

xcp  (semichords) 

0 

4.85 

-0.531 

0.1 

4.84 

-0.531 

0.2 

4.80 

-0.533 

0.3 

4.71 

-0. 535 

0.4 

4.60 

-0. 536 

0.5 

4.44 

-0.542 

0.6 

4.21 

-0.550 

0.7 

3.86 

-0. 563 

00 

• 

o 

3.  36 

-0.578 

0.9 

2.  54 

-0.605 

1.0 


0 


-0.646 


TABLE  4.  SUMMARY  OF  RESULTS 
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Figure  1.  V-g  Curves  For  Case  1  (  Reissner-  Stevens  Theory ) 
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Figure  2.  V-g  Curves  For  Case  2 


r 


Structural  Damping  Coefficient, 


Figure  4a.  V-g  Curves  For  Case  4a 
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Figure  4b.  V-g  Curves  For  Case  4b 
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Figure  5b.  V-g  Curves  For  Case  5b 
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